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Abstract. We consider the behaviour of a dielectric fluid-fluid interface in the 
presence of a strong electric field from a point charge and line charge, respectively, 
both statically and, in the latter case, dynamically. The fluid surface is elevated 
above its undisturbed level until balance is reached between the electromagnetic 
lifting force, gravity and surface tension. We derive ordinary differential equations 
for the shape of the fluid- fluid interface which are solved numerically with standard 
means, demonstrating how the elevation depends on field strength and surface 
tension coefficient. In the dynamic case of a moving line charge, the surface of an 
inviscid liquid-liquid interface is left to oscillate behind the moving charge after it 
has been lifted against the force of gravity. We show how the wavelength of the 
'"O ' oscillations depends on the relative strength of the forces of gravity and inertia, 

^3 ' whereas the amplitude of the oscillations is a nontrivial function of the velocity 

Q at which the line charge moves. 

u ' 

cn ■ PACS numbers: 01.40.Ha, 47.65.-d, 41.20.-q, 47.85.Dh 
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1. Introduction 

(N 

^^ ' Electromagnetic manipulation of dielectric soft matter surfaces has emerged to be a 

r~^ , topic of considerable interest in later years, both from fundamental and applied points 

f^ ' of view. The integration of photonics and microfluidics, optofluidics, is a new analytic 

field - cf., for instance, the reviews [1~3 . 

One of the first examples of this kind was the classic experiment of Ashkin and 

Dziedzic ^ on the elevation of a water surface when illuminated vertically from above 

by a pulsed laser. The force per unit area acting on a dielectric boundary in the 

presence of an electric field is 



X 



Pe = -f ^'Ve, (1) 

The surface elevation of the surface was found to be small, of the order of 1 //m, due 
to the strong surface tension between air and water [i] . Another noteworthy example 
on the interaction between light and matter is the experiment of Zhang and Chang [S] , 
showing how the surface of a micrometer water droplet when illuminated by a laser 
becomes distorted, especially at the rear end because of the lens effect of the sphere. 
A newer series of experiments is provided by the Bordeaux group |B]-[9 , working with 
a two-fluid system of surfactant-coated nanodroplets in oil microemulsions near the 
critical point. The surface tension a can then be reduced to less than a millionth of 
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the usual air-water tension, and the surface displacements can accordingly be very 
large, almost 100 fj,m. The classical electromagnetic theory of this sort of experiments 
is given in several papers; cf., e.g., [5] and the review article |10) . 

The above remarks are related to optical manipulation of liquid surfaces. A close 
sibling of optofluidics is to make use of electrostatic fields instead of optical fields. This 
field of research, often termed electrohydrodynamics, is not new — cf., for instance, 
Taylor's seminal 1969 paper [11] — yet in view of the considerable attention given 
to the optical techniques it is still of interest, not least pedagogically, to understand 
the electrostatic case. Naturally, the electric force density acting in a boundary layer 
is basically the same as above, viz. Eq. ([1]). Since an optical period is much shorter 
than the mechanical response time of a fluid, optofluidic theory works with time- 
averaged fields (over an optical period), and similarities are many. But there are 
obvious differences also. For one thing, the permittivity of a dielectric fluid is usually 
strongly frequency dependent, and is often much higher in the static case than in the 
high-frequency optical case. 

Formation of liquid columns and extraction of jets by electric fields were studied 
in seminal experiments of Raco [12] and Taylor [11], and have been a source of 
investigation at intervals since then. Electrically supported liquid columns have been 
subjected to stability investigations both experimentally and theoretically |13fll6) . and 
the formation of such columns rising against the gravitational force have been shown 
experimentally [iTllTS] . 

The experiment of Raco [12j . in particular, is extreme in that it is very nonlinear 
and moreover is influenced by the finite conductivity of the liquid. In order to describe 
such an experiment, a heavy computational machinery is likely to be necessary. 
However, a useful insight into the fluid-field interaction may be obtained also by 
analytical methods, on the basis of quite simple models. Our intention is just to 
describe the fluid interface response with respect to a strong electric field by taking 
an external point charge Q, or a charged line with charge A per unit length, to be the 
source. This simple model has a high degree of flexibility and its analysis is physically 
transparent as well as mathematically tractable. To some extent, such a theory is 
related to that given by de Gennes [12] on the pushing of a micro droplet using a laser 
beam. 

We will assume in the following that the typical height h of the surface elevation 
is much smaller than the distance L between the charge and the undisturbed surface. 
The reason is that, while the equation giving the fluid deformation in a given electric 
field is simple to derive, the complementary problem of how the electric field from a 
point charge is altered by the presence of an arbitrarily shaped dielectric surface is a 
more difficult problem best solved within computational physics. As is so often the 
case, however, the problem is analytically tractable when restricted to leading order 
in the small quantity h/L, which allows calculation of the field at an arbitrary point 
by method of image charges as if the surface were at rest. Inspection of the governing 
equations which we derive in the following sections is enough to see that including 
corrections to the electric fields due to the surface perturbations provides a higher order 
correction. For consistency we must also linearize the mechanical equations of motion 
for the surface. The problem is nevertheless nontrivial and physically interesting, and 
illuminates key features in the theory of electric liquid surface manipulation. 

Note finally that it is safe to ignore non-linear dielectric response due to the 
quadratic electro-optic effect in the following. Cf. e.g. discussion on p. 147 of Ref. |20) . 
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2. Surface elevation arising from a point charge 



Q (Q") 
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Figure 1. The geometry considered in section [2] 

Consider the setup sketched in Fig. 1: a point charge Q is situated a distance L 
above the undisturbed liquid surface (the xy plane). In cylindrical coordinates, the 
position oi Q is r = 0, z = L. The permittivity in the upper region is ei, and in the 
lower region £2. We assume that 62 > £1 so that the surface force acts upwards. Now 
introduce the image charge 

^, £2 — £1 



-Q, 



(2) 



£2 + El 

situated at r = 0, z = —L. The potential 0i(r, z) at an arbitrary field point P in the 
upper region is that of the charge Q plus the fictitious charge Q' in a uniform medium 
of permittivity £1 [21j . 

where R and R' are the distances QP and Q'P, 



R = ^r^ + {L-z)2, R' = V^2 + (i + z)2. (4) 

Furthermore, if Q" is the image 

Q" = ^^Q, (5) 

situated at the same position as Q, the potential 02 in the lower region is that of Q" 
in a uniform medium of permittivity £2, 

6 ^J-^ (6) 

^"^ 47r£o£2i?' 

These expressions satisfy the boundary conditions at z = 0, namely that (p as well as 
of £d(j}/dz are continuous. 

The radial component Eir = —d(j)i/dr and the vertical component Eiz = 
—d4>i/dz of the field in the upper region are 



47r£o£i V^'"^ -R'V ' 



El, 



1 



47r£o£l 



-Q- 



L-z 



i?3 



Q 



, L + z 



(8) 
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which gives for the square Ef ~ Ef^. + Ef^ at an arbitrary position r, z the expression 



Ef 



1 



(47r£o£i)2 



y , y 2,cjQ 2 _ 1 _ 2\ 



(9) 



Consider now the elevated surface, h — h(r). The electric pressure acting on it is 






— ^In + ^It 
£2 



n 



t/'H'-tJ'" 



n, 



(10) 



where the normal field Ei„ and the tangential field Eit are evaluated on the outer 
surface, z = h~^ , n being the outward normal vector 



n 



-h'r + z f + h'z 



VT+h^ VT+h^ 

Here t is a tangential unit vector and h' = h'{r) is the slope of the surface. Thus 

1 



(11) 



E^t — 



.{Eir+EiM), 



(12) 



VI + /i'2 

which is easily calculable from Eqs. ([7]) and ([U inserting z = h. 

For computational purposes it is convenient to introduce here a reference field 
strength Eq, equal to the field due to Q only, at the origin r = 0, z = 0, 

1 Q 



4:TTeo£i L^ 
We can then write, at position z — h'^ , 

1 , {Q'/Q? , 2Q7Q 



Et = E, 



h^ 



El 



{R/LY {R'lLY {R/LY{R'/LY \ L^ 

(r/L) - (1 - h/L)h' Q' {r/L) + (1 + h/L)h' 



{R'/Lf 



7?2 

" 1 + /l'2 [ [R/Lf 

(observe that El^ ^ Ef ~ Ef^). 

For numerical purposes we use the dimensionless quantities 

r h ^ r - R ~, R' -^ Q' £2 - £i 

L L L L Q £2 + ei 

Noting h' = Ah/Af = h' we may write to leading order in h, h' 



El 



El,^ 



El 



[P + 1)2 

El 



l + Q^ + 2Q- 



~f2_i\ 4(1 -0)~ 



f2 + 1 



r\l + QY + 2f{l~Q' 



(f2 + 1)3 



3r 



-h~h' 



(13) 

(14) 
(15) 

(16) 

(17) 
(18) 



(f2 + 1)3 

plus corrections of order h? ^ hh' . 

The excess pressure force due to surface-tension is related to the two radii of 
curvature i?i and R2 at a given point by |22) 

Ps.t. =Ps.t.n = a(^i- + -i-^n. (19) 

We use the sign convention of '22^ that the curvature is positive when the surface is 
concave upwards, which implies that near the maximum of the elevation, the radii are 
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negative, and the surface tension force points downwards. With a bit of differential 
geometry this may be written in terms of h and its derivatives as [23] 

a d rh' cr , , ,,. , , 

ft.t. = ---=== = - /i' + r/i"+... 20 

r dr Vl + h'^ r ' ' 

where a is the surface-pressure coefficient. The force of gravity is accounted for by 
noting that a surface elevated by a distance h has a hydrostatic pressure which is lower 
than that at the unperturbed surface by a difference 

Aphydrostat. ^ ~{p2 - Pl)gh (2f) 

giving a negative upward (i.e., downward) surface force. Now note that all surface 
forces are directed normally to the surface. The balance equation is therefore scalar 
and the stationary situation is found when the net surface force density is zero. 
Altogether, the governing equation for the surface elevation takes the form 



(P2 - Pi)9h - Ps.t. = 2^o(^2 - £i) 



l-'H'-S'"" 



(22) 



As the density p2 of the lower liquid is larger than the density pi of the upper one, 
both terms on the left hand side of this equation are positive near f = 0, thus both 
adding to the hydrostatic pressure at the level z — 0. This is counterbalanced by the 
field pressure on the right hand side, acting upwards. 
The boundary conditions to equation (j24l) are 



h{oo) ^ h'{0) = 0. (23) 

2.1. Numerical and large-r asymptotic solutions 

Explicitly, thus, the dimensionless surface elevation h is given as the solution to the 
second order differential equation 

h" + [r-' + Afi{r)]h' + [Afo{f) - G]h + As{r) = (24) 



where 



2f( f-£)(l- '^2^ 

(f2 + 1)3 



^i('^) = 7^T—T\3 ' (25) 



J0[r) (f2 + l)4 + (f2+l)3 ' ^^'^) 

(l + Q)V-^+g(l-Q)^ 

^W = pTT)^ • (27) 

and dimensionless parameters 

A^'AEl^IlML, G=^; £-=^. (28) 

2(7 l^ 62 

As is typical for problems involving deformation of liquid surfaces, we may introduce 
the capillary length 



^^^\1 T- (29) 

V {P2- pi)g 

The quantity A is a relative measure of the influence of the electromagnetic force 
compared to that from surface tension. G, similarly, is a measure of the force of gravity 
relative to the surface tension force. In macroscopic geometries surface tension can 
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3 r/L 



Figure 2. Surface elevation for different point charge field strengths (bottom to 
top) Eq = 5 ■ lO^V/m to 1.1 • lO^V/m in steps of lO^V/m. This means Q ranges 
from 5.6 to 12.2 nC. Parameters are L = 1cm, ei = l,e2 = 15, {p2 — Pl)g = 
9790N/m^ . The inset figure is to illustrate the geometry and is not to scale. 



often play a very minor role, so both A and G may take values much greater than 
unity. For a water-air interface, Ic ~ 2.8nim 22J. 

We now turn to the behaviour of h(r) for large r. Letting f — )• oo we have 



l + Ah^l, Afo 
r r 



A(l + Q)2 



G^-G, Asr^ "'^ J. '^' , (30) 

so if we assume h falls off as a power oi r, h oa f^^, then fih' oc h" ex f^^^^ <C foh, 
so these terms can be neglected in Eq. ([M]) . The asymptotic solution is then trivially 
found as 



/i(f) 



A(l + Q)2 



£q£iEq 



£2 - El 1 



f:i' 



(31) 



Gi"^ g{p2 - pi)L 82 + ei r-* 

The asymptotic solution is independent of a as one would expect in an area where the 
elevation varies slowly with r. It is determined by the ratio A/G which is the relative 
influence of electromagnetic forces and gravity, which is also physically obvious. For 
numerical purposes it is more useful to use this asymptotic expansion as a boundary 
condition than h{oo) — 0, since the latter requires needlessly large values of r be used. 
The numerical solution of equation (|24p with boundary conditions (j23p is readily 
found, and shown in figure [21 We use straightforward Runge-Kutta integration with 
a shooting loop to find the value of h{0) so that the asymptotic value (|3T|) is met at 
some large value of r/L. 



3. Surface elevation arising from a line charge 

Consider now the analogous situation where an infinite charged straight line lies 
parallel to the undisturbed surface (the xy plane). The line is parallel to the y axis, 
and intersects the zx plane in the point x — 0, z — L. We let now A denote the charge 
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per unit length. The images A' and A" are related to A exactly as Q' and Q" were to 
Q in the previous section, by expressions © and ([S]). The potentials are [H] 

/>i=-T^ hi-R-- Ini?', (32) 



27reoei 
A' 



27reo£i 



hiR, 



(33) 



27reoe2 
where now R and R' are the distances in the xz plane, 

R = y/x^ + {L~z)\ R' = y/x^ + {L + zy. (34) 

We let henceforth r] — rj{x) denote the surface elevation. The electric pressure P^; is 



given by the same formula (|10|) as before, where now E^ is the tangential component 
of El at the upper surface in the zx plane. An analogous calculation to the one above 
gives for the desired field quantities, at z = ry^. 



E^ ^E^ 



1 



A'' 



1 



2A7A 



(i?/L)2 A2 (R'/L)-^ {R/Lf{R'/LY 



rj- 



L2 



-1 



E- 



El 



(l + ?7'2) 



{x/L) - (1 - ri/L)7f , A' {x/L) + (1 + i^/L)r,' 



{R/Lf A {R'/LY 

where 77' = drj/dx and we have introduced the reference field strength 

Eq = 



27reoei L' 
in analogy to Eq. (fT3|) . 

The surface-tension pressure is now 



Ps.t.n = a- 



V 



n ■ 



{<Jv" + -)n, 



(35) 
(36) 

(37) 
(38) 



(l + r/'2)3/2 

thus a simpler expression than the previous expression (|20l) . The governing equation 
thus becomes very similar to the one we found for the point charge 



(P2 - PljgV - Ps.t. = 2^0(^2 - El) 



£2 



£2 



Ei 



(39) 



with Eqs. ([55|) and ([55]) to be inserted on the right hand side. Boundary conditions 
are as before, 

ry(oo) = r7'(0) = 0. (40) 



3.1. Numerical and asymptotic solution 

For numerical purposes we again introduce the non-dimensional quantities (|16p and 
(PSj) . After linearising equation ([5^ it now reads 

f," + Ah{x)f,' + [AU{x) - G]fi + As{x) - (41) 

with fj = rf/ L, X — x/L and 

2(1-£)(1-A2)i 



Mi) 
six) 



(£2 + 1)2 

2(l-A2)(2x2-F£-a2) 

(x2 + 1)3 ' 

(l + A)25i2+g(l„A)2 
(x2 + 1)2 



(42) 
(43) 
(44) 
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Figure 3. Surface elevation for different line charge field strengths for the same 
parameters as in Fig. [2] and the same values of Eq. This corresponds to A now 
ranging from 5.6 to 12.2 nC/cm. 



Here A = A'/A = Q. 

For large x we proceed as before. We have as x — > oo 

As ^ ^2 ; Afo-Cr^-C; Afi^ Tn (45) 

X^ X"^ 

so, assuming fj oc x^^, fiTJ' <^ ry" <^ ij, we again have a simple asymptotic expression 

A(l + A)2 _ EqEiE^ £2 - El 1 



V' 



-3T, C^OO. 



(46) 



Gi2 g(p2 -/3i)ie2 +ei i^' 

Example numerical calculations are shown in figure |3] for different electric field 
strenghts, calculated in exactly the same way as for the point charge case. 



3.2. The role of surface tension 

Before we move on to the dynamical case, let us briefly discuss the role of surface 
tension. The effect is well illustrated by figure 0] where we plot the elevation with a 
line charge with four different values of the capillary length: once, twice and thrice that 
of water as well as no surface tension. The surface tension force seeks to smoothen the 
surface, making the local radius of curvature as large as possible. Near the maximum 
of the elevation curve, surface tension seeks to pull the surface down, whereas the 
direction of the surface tension force is opposite in the regions where the curve is 
concave up. In general, surface tension counteract any tendency of a surface to 
curve. This is the reason why only small deformations were observed in experiment [4] 
where surface tension was important, as opposed to the laser experiments of Delville's 
group [BHS] who are able to decrease surface tension by many orders of magnitude by 
operating close to a critical point and thereby effectuate large and strongly nonlinear 
interface deformations. 
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4 x/L 



Figure 4. Surface elevation for different values of surface tension. Together with 
the zero surface tension case we plot the results of assuming once twice and thrice 
the capillary length of water, ic, water = 2.8mm. Parameters are otherwise as in 
figure|3with Eq = 10«*V/m. 



4. Dynamical response in the two-dimensional case 

The dynamical response of a hquid-hquid interface to a time- varying radiation pressure 
is of obvious physical interest. There are several papers dealing with this topic, 
theoretical and experimental, when the laser illumination is fixed in space. In addition 
to the references given in Section 1, the recent, extensive treatment of Wunenburger 
et al. [21] is a key reference in that context. 

A natural extension of the kind of static theory laid out in the previous sections 
is to consider the case where the exterior source is not fixed in space, but is moving 
with a constant velocity V in the x direction, parallel to the interface (the xy plane) 
at a vertical separation L. As mentioned above, such a situation has already been 
considered by de Gennes [19] for practical, droplet-pushing purposes. As it is easiest 
to envisage the situation both conceptually and mathematically in the two-dimensional 
case, we will assume that a charged line, directed in the y direction, is the moving 
source. That is, we adopt the same physical model as in Section [31 Our purpose is to 
examine the time-dependent elevation of the interface, for a given value of V. We will 
for simplicity assume zero viscosity in the liquids, which means that there is no wake 
formation behind the line moving through medium 1. We assume V is non-relativistic, 
and that the magnetic field from the moving charge is thus negligible compared to the 
electric. 

In the discussions in this chapter we shall omit the effects of both viscosity 
and surface tension in order to make the physical picture more transparent and also 
simplify the mathematics somewhat. Inclusion of both these physical corrections is 
straightforward using the same analytical techniques as presented herein, however, 
and could form a useful extension of it. 

Before we commence, we make the physically natural ansatz that all physical 
quantities depend on x and t only through the Galilei transformation variable 
^ = X — Vt. This is the same as assuming that in a coordinate system where the 
charge is at rest, the surface elevation appears stationary. With this we may replace 
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cv 




a+ 


'i = 




a- 











Figure 5. The geometry of the small control volume used to derive the equation 
of motion in the dynamic case. 



-w 



(47) 



in the following 

d__ d_ d_ 
9^~ df' di 

Consider now a tall and thin control volume CV as shown in figure [5l intersecting 
the interface with a small end area A whose extension is so small that ri{(^) is 
approximately linear inside CV. We let 77 mean the elevation in the center of CV. 
Newton's second law for such a volume, which follows from Reynolds' transport 
theorem, asserts that the z component of the sum of external forces (here pressure 
and gravity) equals the change of fluid momentum within CV (e.g.. Chapter 3 of [25'): 



E 



-4 



(48) 



F^""' = f pwv ■dA+— pwdV 

Here we let the decomposition of the velocity vector be v = ux + wz, dA is directed 
out of the control surface 9CV, and dV denotes volume integration. 

The sum of the external forces, including gravity, hydrostatic pressure and the 
pressure discontinuity at the surface due to electromagnetic forces and surface tension, 
is easily ascertained to be given by the same expressions as used in the previous 
sections, 

Y, FT' = -4(PEM + Ps.t. - ^P9V) (49) 

and we need not derive this again, pem and Ps.t. have the same meaning as in section 
[3J but with ^ replacing x in pem . 

Treating the right hand side of equation P5|) is a little more involved. First 



let 

the heights a~ and a+ (see figure [S]) be large, so that w is zero at the end surfaces 
(clearly, w must tend to zero far from the interface). To linear order there is no flux 
across the side surfaces of dCY. That disposes of the surface integral. We are left with 



■Fr{o^Ap2 



w{^, z)dz + Ap_ 



d_ 
di 



w{^,z)dz. (50) 



We now require knowledge of how the z- velocity w depends on x and z. Since the 
surface elevation is assumed small, we are free to make use of the velocity potential 
(j), where v = V0. Assuming incompressibility gives the continuity condition 'V'^cj) = 
which upon further differentiation w.r.t. z implies 

-jU^O. (51) 



\dx'^ 



dz 
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As boundary conditions we clearly must demand that w{^, z = ±c») — 0. Furthermore 
the so-called kinematic boundary condition at a surface is 

d77 

w{^, z = T]) = fj + U-— w ?7, (52) 

dec 

since the second term in the middle form is of quadratic order in rj/ L. 

According to the standard method of separation of variables we assume 'w{£,, z) — 

5(^)Z(z) and conclude that S" = —k'^E and Z" — k^Z for some real number k. In 

Fourier form the solution may then be written 

w{^,z) = VL u;(fc)e-'=l^-"le'^-« (53) 

Jo 27r 



which is seen to solve (|5ip and fulfil boundary conditions at infinite |z|. Of course, 
the physical value is the real part of this. We next define the complex elevation 

VciO^L'J^ ^V{k)e"'^ (54) 

so the kinematic boundary condition ([5^ gives w{k) = —\kLrj{k). With the form (15^ 
the integrals in equation (j50p are easy to solve, and inserting the answer into the left 
hand side of (|49l) after differentiating using the ansatz (|47| we obtain the equation 

PEM(C,f7c,?7c) +Ps.t.('7c) - ApOTc 
r°° Ah 
^-V^L\p,+p,) fc^(fc)e"=«=iF2(pi+p2)r?^ (55) 

Now we may let rjc = rj + irji and write down the real and imaginary parts 
of equation ([55]) separately. To simplify the picture somewhat, consider the system 
when there is no surface tension, a — 0. For the parameters used in the previous 
sections, this hardly changes the results, but makes the physics more transparent and 
also simplifies the numerical treatment considerably since stiffness is avoided and the 
order of the differential problem is reduced from 4 to 2. Eliminating rji from the 
equations we obtain our governing equation without viscosity or surface tension, 

(/i^ + /o + r) (/i,?' + fov + s - Fry) - -T^". (56) 

In order to have nondimensional parameters not containing a, we have defined here 

G ^ 2gL{p2 - pi) . ^^ 2V^{p2+Pi) 
A eo{e2-ei)E^' eo{e2-ei)E^' 

Obviously these parameters are measures of the relative strength of gravitational to 
electromagnetic forces, and inertial to electromagnetic forces, respectively. Explicitly, 

if! + r')fl" + ififi + 2/1/0 - 2F/i)77' 

+ [/i/o + (/o - r)2]ry + f,s' + (/o - F),s - 0. (58) 

where /i(Oi /o(0: ^(0 (argument S, suppressed above) are as defined in (li^ - pi)) upon 
X — )■ ^. We have nondimensionalized as before ^ = ^/L. 



'^ - A - ^ f^ , NZ72' '^ - ^^(^^_^\J^2- y^') 
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4.1. Asymptotics and numerical solution 

Repeating our previous exercise of taking the large ^ asymptotics of the different 
coefficient functions we find that the asymptotic form of 77(^) is unchanged from the 
static case if we assume the solution to fall off as a power law as ^"'^ at infinity, 

r g[P2- pi)Le2 + £i^^ 

The asymptotic form of fj' is thus 77' ^ -~2fj/£^^, which together form a useful set 
of boundary conditions at large £,. This condition must clearly hold in front of the 
moving charge, in accordance with our assumption that the liquid surface is at rest 
before the arrival of the line charge. 

However, since our liquid is assumed inviscid, fj no longer needs to tend to zero 
at infinity now that we consider a system which is not in mechanical equilibrium. To 
wit, behind the moving line charge we might expect waves to form since the charge, 
after setting the surface in motion, moves on and leaves it to oscillate in its wake. 
When not assuming ry to tend to zero, the asymptotic equation from ()58|) must also 
include a 77" term, 

TV + r^^^-r^ii^. (60) 

We shall not dwell on the particular solution to this equation (which can be found!), 
but note rather that equation (1601) has two homogeneous solution which do not vanish 
at infinity, 

77 = Ci cos — + C2 sin — (61) 

with integration constants C'l and C2 which we are not in a position to evaluate in a 
simple way. This asymptotic form must be valid far behind the moving charge where 
the surface has been disturbed and moves as an undamped wave. 

Notice at this point how our theory is valid when V = (the static case considered 
before) and when V is sufficiently large, but not when V/^/gL is small. The reason is 
that we have assumed not only that ij -^ 1, but also 77' ^ 1, in order to get rid of the 
terms (1 + 77'^)^^/^. However, when V is small, the wavelength of the waves in the 
wake becomes short, 

27TT V' P2+PI .„„. 

A = -— =27r— 62 

r gL p2- pi 

and the derivative becomes proportional to gLV^^. Note that this wavelength is 
independent of the electromagnetic properties, and is simply a measure of the relative 
strength of the forces of inertia and gravity. 

In figure [6] we plot the solution for a line charge moving towards the right at 
different velocities. At large V there are no problems, but already at the two lowest 
values of V we are in fact pushing the boundaries of linear theory, since 77' goes as 
high as 0.5. Equation (|58p predicts very high maximum elevation for small velocities, 
but this is wrong and happens because the compensating factor, such as the square 
root denominator in equation (1551) , is not included. The highest peaks may be slightly 
overestimated in figure [51 therefore. 

We observe a number of interesting features. For example, when velocity is small, 
the peak is in fact higher than the static, mechanical equilibrium value. This is because 
the rising liquid continues beyond the equilibrium value before the influence of gravity 
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Figure 6. Surface elevation in the absence of surface tension and viscosity for a 
charge moving at velocity V towards the right [particle position {x, z) = {Vt, L)\. 
We let p2g = 9790N/m^ and pi fa as in the case of a water/air interface. The 
dotted line shows the static mechanical equilibrium situation (also for ct = 0). 



finally makes it turn and fall back again. On the other hand, when the line charge 
moves fast, the fluid only reaches its maximum velocity some time after the charge 
has moved on. We also see that the amplitude of the harmonic oscillations is a non- 
trivial function of V . When the movement is very slow, oscillations are very small, 
since the charge, after lifting the surface, stays around to set it gently back down 
again. This is what we see in figure |6] for the two slowest velocities. When V is very 
high, likewise, oscillations must, in addition to having long wavelength, be of small 
amplitude, since the surface hardly has time to feel the presence of the charge before 
it is gone again. Somewhere in between these extremes is found the velocity which 
gives the maximum oscillation amplitude. Finally, we see that for small velocities, the 
first peak appears in front of the line charge, as is the case for V = {)A\/gL in figurelHl 
For even smaller velocities than this, the surface is lifted up, then drops again due to 
gravity, whereupon it is lifted again to a second time, since the field from the charge 
is still present and strongly felt. The beginnings of a second peak are observed in 
the OA^/gL graph (to fully study this phenomenon, nonlinear theory must be invoked 
as previously discussed). All in all, our simple set-up treated with linear theory has 
yielded surprisingly rich physics. 

In the previous we have omitted the effects of viscosity, although it is obvious 
what the effect of its inclusion would be. We have seen that far behind the moving 
charge, the governing equation takes the form of an undamped harmonic oscillator (or 
wave) equation. A viscosity term, being proportional to fj would include a damping 
in this equation, making sure the harmonic waves be damped exponentially behind 
the moving charge, admitting overdamped and underdamped solutions depending 
on the relative value of viscosity and wave number F/T. Studying this effect is 
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straightforward within the above framework [a friction force is simply added to the 
right hand side of equation (|49p ]. Along with the effects of surface tension it serves as 
a useful extension problem for pedagogical purposes. 

5. Conclusions 

We have considered the basic electrohydrodynamic case of a point charge and a line 
charge above a fluid-fluid interface. The electric field from the charge acts on the 
dielectric interface, causing it to rise above the equilibrium level. In the static case 
the shape of the fluid-fluid interface is determined by the relative strength of three 
forces acting on it: the electromagnetic lifting force, gravity and surface tension. 

For small elevations the deformation can be described by a single ordinary 
differential equation which is readily solved numerically. We derive and solve this 
equation for the case of a point charge and a line charge. When surface tension is 
small the maximum height of the surface is determined by a dimensionless number 
which is the relative weight of the electric lifting force vs. gravity. The effect of surface 
tension is to counteract the tendency of the other forces to induce surface curvature, 
which we illustrate in figure |4l 

Finally we derived the dynamical equation of motion for a surface subject to a line 
charge moving at a constant velocity. Altough apparently just a small generalization 
of the static case, the system is surprisingly rich in physics. As the moving line 
charge goes past, the liquid surfaces is elevated from its unperturbed level and is 
left to oscillate around an equilibrium level as the charge moves away [we assume 
inviscid fluid, so the oscillations are undampled] . The details of the time dependent 
movement of the fluid surface, however, depend nontrivially on the speed at which the 
charge moves, compared to a reference velocity y/gL [L being the smallest distance 
from charge to undisturbed surface] . For high velocities the surface hardly has time to 
react to the presence of the charge, and the oscillations have long wavelength and small 
amplitude. For low velocities the oscillation wavelength is short, and the maximum 
elevation can be considerably higher than that of static equilibrium. The amplitude of 
oscillations behind the charge are small also in this extreme, however, since the charge 
has time to set the surface gently back behind it. In between these extremes is found 
the velocity giving the highest amplitude of oscillations behind the moving charge. 
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